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ABSTRACT 

We present an analysis of two thermodynamic techniques for determining equilibria of self- 
gravitating systems. One is the Lyndcn-Bell entropy maximization analysis that introduced violent 
relaxation. Since we do not use the Stirling approximation which is invalid at small occupation num- 
bers, our systems have finite mass, unlike Lyndon-Bell's isothermal spheres. (Instead of Stirling, 
we utilize a very accurate smooth approximation for In a;!.) The second analysis extends entropy 
production cxtremization to self-gravitating systems, also without the use of the Stirling approxima- 
tion. In addition to the Lynden-Bcll (LB) statistical family characterized by the exclusion principle 
in phase-space, and designed to treat collisionless systems, we also apply the two approaches to the 
Maxwell-Boltzmann (MB) families, which have no exclusion principle and hence represent collisional 
systems. We implicitly assume that all of the phase-space is equally accessible. We derive entropy 
production expressions for both families, and give the extremum conditions for entropy production. 
Surprisingly, our analysis indicates that extremizing entropy production rate results in systems that 
have maximum entropy, in both LB and MB statistics. In other words, both thermodynamic ap- 
proaches lead to the same equilibrium structures. 

Subject headings: galaxies: structure — galaxies:kincmatics and dynamics 



1. INTRODUCTION 
1.1. Motivation 

Understanding how collisionless systems attain a spe- 
cific mechanical equilibrium state is fundamentally im- 
portant to astrophysics. For example, the cold dark mat- 
ter structures that exist around galaxies are expected to 
fall into this class of system. Individual dark matter 
constituents (whatever they may be) should evolve ac- 
cording to a mean-field gravitational potential, free of 
the influence of individual encounters. 

The range of mechanical equilibria available to a colli- 
sionless system is defined by the Jeans equation, which 
represents the condition that no portion of the system ex- 
periences a net force. Unfortunately, the Jeans equation 
admits an infinity of solutions. Even if the mass distri- 
bution is specified, there is an infinite set of acceptable 
mechanical equilibria, each involving a different velocity 
distribution. For spherical systems, these velocity distri- 
butions differ in their anisotropy profile that quantifies 
radial versus tangential motion. However, the question 
remains, how and/or why does any one collisionless sys- 
tem evolve to its particular mechanical equilibrium end- 
state, and what are the properties of such a state? 

1.2. Thermodynamic Approaches to the Problem 

The statistical mechanics description of thermodynam- 
ics provides one path to obtaining the description of the 
final relaxed state. A fully relaxed system is the most sta- 
tistically likely state of that system, or the one with an 
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entropy maximum. In calculating the most likely state 
it is implicitly assumed that all states of the system are 
equally accessible. 

Such an approach was taken bv lLynden-Bell (|1967D and 
applied to self-gravitating collisionless systems, with the 
hope of explaining the observed light distribution of el- 
liptical galaxies. A collisionless system can be thought of 
as a fluid in the 6D phase-space of position and velocity. 
The 'particles' in Lynden-Bell's analysis are parcels of 
phase-space, i.e., parcels of this fluid, and so the distribu- 
tion function (DF) representing the phase-space density 
is defined in terms of energy per unit mass, not energy 
per particle. The analysis resulted in a DF similar to 
the Fermi-Dirac ca se, but with a different normalization. 
iLvnden-Belll (|1967f ) argued that a non-degenerate limit is 
appropriate for stellar systems, and thus arrived at a DF 
similar to that of the Maxwell-Boltzmann case (an ex- 
ponential) which resulted in the isothermal sphere repre- 
senting thermal equilibrium. Since the isothermal sphere 
has an infinite extent and mass, its emergence from the 
entropy maximization procedure, which demanded a fi- 
nite mass system, presented a contradiction. This appar- 
ent failure of entropy maximization was puzzling, and it 
was often argued that such systems do not have states of 
maximum entropy. Some effort was made to investigate 
maximizing entropy with additional c onstraints beyond 
mass , energy, and angular mom entum (|Stiavelli k Berth] 
119871: iWhite k Naravanlll987t) . Other routes involving 
minimum energy sta tes of self -gravitating systems were 
also develop ed (e.g., IAly||1994D . 

Recently. | Madsenl (| 1996T) (based on earlier work by 
I Simons! 119941) has pointed out that the reason for the 
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system's infinite mass was the use of the Stirling approx- 
imation, Inn,! w nlnn — n. In contrast to systems usually 
treated in standard statistical mechanics, self-gravitating 
systems can apparently have small phase-space occupa- 
tion numbers n, making Stirling a poor approximation. 
Specifically, these systems have regions of phase-space 
or energy-space that are nearly or completely unoccu- 
pied, such that n will be small. Spatially, these regions 
can correspond to the center of the potential as well as 
its outer edge. Using Maxwell -Boltzmann statistics and 
the exact Inn!, iMadsenl (|1996[ ) has found a distribution 
fu nction from entropy maximization that is very s imilar 
to iKind (|1966l) models. iHiorth fc Williams! (|2010D have 
shown that the energy-space occupation function N(E) 
derived using a very accurate smooth approximation to 
Inn! closely resembles the results of collisionless TV-body 
simulations (jWilliams et aLll2010( ). 

1.3. Brief Review of Statistical Representations 

From a statistical point of view, entropy is simply a 
measurement of the number of states accessible to a par- 
ticular system. This relationship is most commonly ex- 
pressed quantitatively as, 

5 = fc B ln^, (1) 

where f2 is the number of accessible states and fee is 
the Boltzmann constant which serves to give entropy 
the correct thermodynamic units. As a result, count- 
ing procedures are key to determin ing specific realiza- 
tions of entropy. iLvnden-Belll (|1967t ) discusses how there 
are four counting types that lead to physically relevant 
situations. Bose-Einstein statistics follow from counting 
states for indistinguishable particles that can co-habitatc 
in the same state. When indistinguishable particles are 
not allowed to share states, Fermi-Dirac (FD) statistics 
emerge. Classical Maxwell-Boltzmann (MB) statistics 
result from counting states available to distinguishable 
particles that can share states. Completing the symme- 
try, systems where distinguishable 'particles' — actually, 
parcels of phase-space — cannot co-occupy states obey 
what has become known as Lynden-Bcll (LB) statistics 
(statistics 'IV in Lynden-Bell's original notation). Each 
type of statistics will produce different representations of 
entropy, but we will focus on the two that deal with clas- 
sical, or distinguishable particles, namely, LB and MB. 

We briefly recap the notation used in ILvnden-Belll 
(|1967t ) before proceeding with our discussion. The six- 
dimensional position- velocity phase-space (x,v) is the 
usual setting for determining the statistics. Imagine 
phase-space to be divided into a very large number of 
nearly infinitesimal parcels, called micro-cells, each hav- 
ing volume w. Each micro-cell can either be occupied 
or unoccupied by one of the N phase-space elements of 
the system. These elements can be thought of as repre- 
senting the fine-grained distribution function, which has 
a constant density value n. Because collisionless pro- 
cesses imply incompressibility of the fine-grained distri- 
bution function, the phase elements cannot co-habitate. 
If phase-space is also partitioned on a coarser level so 
that some number v of micro-cells occupy a macro-cell, 
then we can discuss a coarse-grained distribution func- 
tion. The volume of a macro-cell is then vvo and the ith 
macro-cell contains m phase elements. We assume that 



while the volume of a macro-cell is much larger than 
that of a micro-cell, it is still very small compared to the 
full extent of phase-space occupied by the system. The 
number of ways of organizing the n»s elements into the v 
micro-cells without co-habitation is, 

If the elements were allowed to multiply occupy micro- 
cells, as in MB statistics, this number would be given by 
v ni . 

To get the total number of accessible states, the pos- 
sible ways to distribute the N phase eleme n ts int o rij 
chunks must also be included. iLynden-B cll (1963) de- 
rives, 

In the MB case, the only change is that the factorial ratio 
in the final product term is replaced by v ni . 

The LB case disallows two phase-space elements from 
inhabiting the same phase-space location, so it explicitly 
takes into account the incompressibility of a collisionless 
fluid. Since we are primarily interested in dark matter 
halos, LB is the natural case to consider. For complete- 
ness, and for the sake of having a comparison, we also 
treat the MB case, in the Appendix. 

We argue that the lack of an exclusion principle in 
the MB case is equivalent to allowing collisions between 
particles. In a collisional system, particles from distant 
phase-space locations can be scattered into any other 
phase-space location, thereby increasing the phase-space 
density at the latter location. In principle, there is no 
limit to how high the density can get through such scat- 
terings. In practice, the phase-space density probably 
can not become very high at most locations, but it can 
be higher than the original fine-grained DF. We note that 
it is common to use MB to r epresent collisional systems. 
For example. IMadsenl (|1996[ ) argues that it is the correct 
statistics to use for globular clusters where the relatively 
small number of stars allows the cluster to relax through 
two-body interactions. It then makes sense that the en- 
ergy distribution that the cluster will arrive at will be the 
same as that in a cloud of gas, which relaxes through col- 
lisions between molecules. In the non-degenerate limit, 
when the micro-cells are very sparsely populated and the 
density of the coarse-grained distribution function is very 
dilute compared to that of the fine-grained function, both 
LB and MB distribution functions, and hence density 
profiles, will look the same. 

1.4. This Work 

In this paper we explore two possible approaches to de- 
riving the final equilibrium state of self-gravitating sys- 
tems, for each of the two types of statistics, LB and MB. 
In both, we use a very accurate, smooth approximation 
for In a;! valid for arbitrary occupation numbers x. How- 
ever, the price we pay for this improved approximation 
is the loss of analytic solutions. 

The first approach assumes that the final state is the 
maximum entropy state, an assumption that was first 
used in the context of self-gravitating systems in 1950's 
(|Ogorodnikovl ll957) . The second approach, again in the 
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context of self-gravitating systems, was first taken by 
iBarnes fe~ Williams ([201 ff ); it posits that the final state 
corresponds to the extremum of entropy production rate. 

We explore cxtrcmizing entropy production because 
it has not been established beyond a doubt that real 
or computer simulated systems do fully relax to max- 
imum entropy states. It is possible that their steady- 
state configurations do not correspond to maximum en- 
tropy. Prior work on thermal non-equilibrium systems, 
but not in astrophysical contexts, suggests that station- 
ary states — like mechanical equilibrium — occur when en- 
tropy production is extremized (e.g., iPrigogiiiel 119611 : 
iJavnesI ll~980l : Ide Groot fe Mazuri fl~98l IGrandvl 120081) . 
We investigate their applicability to self-gravitating sys- 
tems in mechanical equilibrium. 

The new aspect in the present paper is that we use a 
very accurate approximation for lnx!, unlike our previ- 
ous pap e r that assumed the Stirling approximation. As 
iMadsenl (|1996f ) has shown, replacing Stirling with an ac- 
curate approximation (i) results in systems with finite 
total mass and energy, and (ii) significantly changes the 
structure of the systems. 

In all, we present four derivations; entropy maximiza- 
tion for the LB a nd MB statistics are covered in Sec- 
tions 13.11 and IA.11 Extrcmization of entropy production 
for L B an d MB statistics are carried out in Sections 13.21 
and IA.21 respectively. We develop expressions for the 
relaxation functions (see § H]), and use these to better 
understand the evolution of coarse-grained distribution 
function. We compare our results with analogous ver- 
sions of entropy producti on derived using th e stan dard 
Stirling approximation in IBarnes fc Williams! ([2011D . 

Figure[T]puts the present paper (BWII in the figure) in 
context. It is a schematic summary of the various statis- 
tical mechanical approaches to self-gravitating systems. 
The possible ways to frame the problem appears at the 
top of the figure; one can formulate the problem in ei- 
ther the regular phase-space, or the energy space. Below 
the thick horizontal line we show the two different routes 
for attaining the final steady-state state: maximizing en- 
tropy, and extremizing entropy production. Once these 
choices are made one has to decide whether small oc- 
cupation number regime will be important or not, and 
hence whether to use the Stirling approximation for In a;!, 
or not. In the latter case, one must then decide whether 
to use the discrete ("discr.") step-like, i.e. exact version 
of lnx!, or to approximate it with some smooth function 
("cont.") which remains very accurate down to small 
x. Note that HW10 and the present paper use different 
but similar approximations. There is no physical reason 
to introduce the exclusion principle in the energy state- 
space, hence the corresponding regions are marked as 
"not relevant" . The bottom entries of some columns in 
the table contain names of papers where the correspond- 
ing options were considered. K66 in parentheses below 
BWII means thatlKing ( 1966;) results are ne arly identical 
to our s. M = oo under LB67 means that iLvnden-Belll 
(|1967t ) final result, the isothermal sphere, had infinite 
mass. 

Much of t he background mat e rial fo r this work may 
be found in IBarnes fcW illiams (|201lD , and we briefly 
summarize these previously obtained results in Section[5J 



SUMMARY OF RESULTS FOR LARGE OCCUPATION 
NUMBERS 



IBarnes fe Williams! (|2011D have investigated entropy 
production in self-gravitating systems described by MB 
and LB and have developed expressions for the entropy 
production, a for both the MB, 
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and LB cases, 



0"LB 



wq 



In 



1 - In N 



dv. 



V - / 



o 



dv, 



(4) 



(5) 



where constant C in Equation [5] is (IniV — 1) + 
(1/^)^ ^111^. In these expressions, / is the coarse- 
grained distribution function and T is the relaxation 
function and forms the right-hand side of the Boltzmann 
equation, 



dt 



Vf + a-V v f = T(f). 



(6) 



If one were to assume, for a moment, that in the above 
equation / is a fine-grained DF, then the right hand side 
would be zero for collisionless systems. In other words 
r = 0, which means that on fine-grained scales there is 
no change in, or production of entropy; in collisionless 
systems, entropy is fixed throughout evolution. Let us 
be clear, we are not advocating for a specific process as 
a source for T, like collisions in a gas. The relaxation 
function simply describes the Lagrangian time rate of 
change of the distribution function. Returning to our 
case where / represents the coarse-grained DF, Equa- 
ti on [HI states that entropy is produced [and as we argue 
in IBarnes fc Williams! (|2011[ ) it happens even in systems 
that have attained macroscopic steady-state] because on 
microscopic scales the fine-grained DF continues to wind 
and twist, which when combined with coarse-graining, 
gives rise to no n-zero entropy c hange. From a different 
starting point, iChavanisI (|1998[ ) develops an expression 
for the right-hand side of Equation [5] in terms of a "dif- 
fusion current" that relates to correlations between fluc- 
tuations in the fine-grained distribution function. While 
that work details the makeup of this diffusion current, 
we simply focus on the broad behavior of the relaxation 
function. 

We find extremum entropy production conditions by 
setting the variation of entropy production, 8a equal to 
zero. This operation gives expressions for the relaxation 
function. For the MB case, 
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(1- In AQr MB (f = y) 
In (//?;) + 1 -In N 



(7) 



The LB relaxation function is slightly more complex, 

-C7r LB (/ = r;/2) 

^ (/) " ln[// („-/)]-C- (8) 

Like lLynden-Bel (fl967h . the|B arnes fc Williamsl (pOlTh 
work assumes that the large n Stirling approximation is 
valid for the systems being investigated. Here, we will be 
deriving relations analogous to Equations HI El El an d [HI 
but using a very accurate approximation, after discussing 
the results of entropy maximization below. 
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3. RESULTS FOR ARBITRARY OCCUPATION NUMBERS 
The approximation that we utilize in this work is, 



In a;! = (x + — ) In (x + 1) 



hl27T 



+ Aq, 



where 



Aq, : 



(z 2 + 2z+f|) 

[X + 12/ 



(9) 



(10) 



In the large x limit, this reduces to the usual Stirling 
approximation In a;! = x In x — x. Our approximation 
is slig htly different from the one in iHiorth fc Williams! 
(|2010D . Comparisons between the Stirling approxi- 
mation, IHiorth fc Williams! (|2010D approximation, and 
Equation [§] are shown in Figure [2 The plots illustrate 
the function ip(x + 1) = d\nx\/dx for the three cases. 
The HW10 and Equation [5] approximations are nearly 
identical. Most importantly, these two approximations 
are well-behaved at x = 0, unlike the Stirling approxi- 
mation. 

3.1. Entropy Maximization 

In this section, we follow the overall path taken in 
iLvnden-Belll JT967) to determine the description of en- 
tropy for the LB statistics family. Briefly, counting ar- 
guments for phase-space macro-cell occupation are com- 
bined with the statistical definition of entropy to give 
specific representations. Variations in entropy assuming 
constant mass and energy are then set to zero in order 
to determine the necessary distribution functions. 

If one demands that multiple phase-space elements 
cannot simultaneously occupy micro-cells, the multiplic- 
ity of states (Equation [3]) combined with the definition 
of entropy results in, 

S'lb = k B In AH _ ^ \ n71i \ + ^ m v \ - ^ m (y - ni )\ 

i i i 

(ii) 

where again the summations run over the number of 
macro-cells. We will assume that both N and v are much 
larger than 1, so that Stirling's approximation is valid for 
use in the first and third terms. However, for the and 
v — Hi terms we will use our improved approximation, 



Inn,! = (m + 1/2) In (m + 1) - m + 



lll27T 



+ Ao, ni . (12) 



Note that this usage does not require either m or v — 
Hi to actually be a small value, rather this keeps the 
accounting accurate in the event that they do become 
small. If these values are always large, there will be no 
difference from the Stirling approximation. The entropy 
expression now reads, 

Slb = S LB ,o -kB^2 [("< + V 2 ) In fa + 1)+ 

i 

{v-m + 1/2) In {y - m + 1)4- X Q , nt + 



A 



(13) 



where S LB ,o = [N In N - N + M(v In v - In 2n)\ , and 
M is the total number of macro-cells. 

Transforming from the discrete macro-cell occupation 
number to the continuous coarse-grained distribution 



function /, the entropy becomes, 



S'- 
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dx dv. 



+ \uf/n + 

(14) 



Taking the variation of this entropy expression to be zero, 
with constant mass and energy constraints, leads to the 
following condition, 



ln(F+ 1) 



OF +1/2) 



(^ + 1) 



-\n(v-F+l)- 



{v-F + 1/2) 
(iv-F+1) 

dAo.F dAo : y- 
dF dF 



+ /i + f3e = 0, 



(15) 



where F — vf ' jr\ is a scaled coarse-grained distribution 
function and [i and j3 are undetermined multipliers asso- 
ciated with mass and energy conservation, respectively. 
The e term is the specific energy of a phase element lo- 
cated at position x with velocity v, e — i; 2 /2 + $. The 
derivative of the A function is, 



dA 0if 
dF 



-(F+l) 
(F + 600/576)' 



(16) 



After substituting for these A derivatives and combining 
terms we have, 



In 



F + l 



f -F + l 



{F-v/2) 



(F + l)(u-F + l) 



2F 2 - 2F^ - (1 + 600/576)i/ - 600/288 

F 2 -Fv- (600/576)i/ + (600/576) 2 + 
H + pe = 0. (17) 

The inelegant ratio terms in this expression are both 
symmetric about F = v/2. The second term on the 
left-hand side of Equation [T7] has values of —1/2 when 
F = 0, when F = v/2, and 1/2 when F = v (the 
maximum value of F for the LB case). The third term 
on the left-hand side of Equation [171 is a nearly constant 
function with a value very close to —2 for < F < v. 

We have not attempted to find an analytic solution 
for F. Graphical solutions of Equation [17] for a series 
of e values produces the picture of f /rj seen in Figure G3 
The overall character of the distribution function is the 
same in the non-Stirling and Stirling versions, and both 
functions are very similar to Fermi-Dirac distribution, 



/) 



FD 



exp — (// + /3e) 
1 + exp -(/! + /3e) ' 



(18) 



This distribution function of Equation[T7]can be trans- 
formed into a density distribution using the Poisson 
equation and the fact that 



p(r) = / f(r,v)dv = 4tt / f(e)y/2[e - $(r)]de. (19) 



Note that this procedure imposes an isotropic velocity 
distribution for the system. 
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The density and logarithmic density slope correspond- 
ing to the non-Stirling function are presented in Figures[4] 
and [5] and will be discussed further in Section |4~T1 



3.2. Entropy Production Extremization 

Just because it is possible to describe thermal equi- 
librium states for collisionlcss self-gravitating systems, it 
does not follow that real systems (either physical or sim- 
ulated) must achieve them in a Hubble time. It is pos- 
sible that real systems incompletely relax, leaving them 
in a t hermal non-equ ilibrium, but long-lived stationary 
state. iMadsenl ()1996l ) also points out that in very slowly 
evolving systems, the maximization of entropy is only 
temporary. This leads us to infer that the production of 
entropy may be a useful quantity for discussing quasi- 
cquilibria of collisionlcss systems. With that in mind, we 
now proceed to develop the conditions required for a ther- 
mal non-equilibrium state to be stationary. Specifically, 
we find an expression for entropy production in the LB 
statistical family and then extremize it. (As for the pre- 
ceding section, an analogous derivation for the Maxwell- 
Boltzmann statistics may be found in Appendix lA.2[ ') As 
discussed in Section 11.41 existing work on thermal non- 
equilibrium in non-astrophysical settings suggests that 
stationary states occur when entropy production is ei- 
ther maximum or minimum (Ide Groot fc Mazurl Il98i 
IGrandvll2008l and references therein). 

From Equation 1141 the entropy density in the Lynden- 
Bcll case can be written as, 



PSLB 



VTU 



(F+l/2)ln(F+l)- 



[y -F + 1/2) In (i/ -F+l 
,Slb,o 



Xq.F + \),u-F — F- 



Nfo 



dv, 



(20) 



where, as before, F = vf/rj. 

Taking a partial time derivative of Equation 1201 results 



9 / x k B f OF n , _ . 
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-\n{v - F +1) 



F + l 

9\o,F , 9Ao, y -F 



v- F+l/2 
v -F + 1 



+ 



dF 



dF 



C 



dv, 



where C = Si^o/A^b ^ s a constant. 
Upon substituting dF/dt from the Boltzmann equa- 



tion into Equation 1211 we get a lengthy expression, 
d 



dt 



{PS-LB) 



H-vVF) 
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In (F+l) 
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(^ + 1) 



ln(,-F + l )+ -- F „ +1 / 2 + ^ + ^-C 
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(-a ■ V V F) 



In (F + l) 



dF OF 

(F + l/2) _ 
{F + l) 



dv — 



]n{v-F+l) 



v- F+l/2 <9A ,i 



d\ 



u-F 



In (F + 1 



v-F + l dF 
(F + l/2) 



OF 



-C 



dv — 



In (y - F + 1) 



(F + l) 
v- F+l/2 d\ ,F 



v-F + l 



dF 



dF 



-C 



dv, 



(22) 



where 7 = vT/rj. 

The first term on the right-hand side of Equation [22] 
can be transformed into —V • ps^v. The acceleration- 
dependent terms in Equation [22] all disappear under the 
assumptions that the distribution function is even in v 
and disappears when v = w max - The final representation 
of the time rate of change of the entropy density is, 

|(^lb) = 

-— V- / [(F + l/2) In (F + l) - 
vzu J 

(v-F+ 1/2) ln(v - F + 1) + Xa, F + X , v -f - FC]v dv 



7 (F) 



In (F + 1) + F+l / 2 _ l n [y _ p + 1) 



F + l 

v- F + l/2 dA ,_F dAo,„-F 



v-F + l 



dF 



dF 



C 



dv. 



(23) 



Again assuming the velocity field to be composed of mean 
and peculiar components v = Vq + v p , we draw a corre- 
spondence between the terms in, 



d_ 

dt 



(PSlb) = - V • pS L B«0- 

V- J[(F +1/2) ln (F + l) - 



(v - F + 1/2) ln(v - F + 1) + A ,f + Xq, v -f - FC]v p dv 



(21) — 



VVD 



In (F + 1) + f + l { 2 - ln (v - F + 1) 



F+l 



v -F + 1/2 dA ,F dAo.^-j 



C 



dv. 



v-F + l dF dF 
and those in the continuous version of the entropy density 



(24) 
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evolu tion equation, (Equation 6 in IBarnes fc Williams! 

EnnD, 

^(ps) = -V ■ + psv ) + a. (25) 

The entropy flux due to random motions S is given 
by the integral in the second term on the right-hand side 
of Equation 1241 The remaining term then makes up the 
entropy production for the system, 



0"LB 



k B 
vw 



ln(F+ 1) 



J 1 +1/2 
F+l 



- In (v - F + 1) 



F+l/2 , dA , F , dAo,„-F 



dF 



dF 



C 



du. 



(26) 



As expected, the relaxation function determines the en- 
tropy production rate for the system. 

We find that the condition for the entropy production 
term in Equation 1261 to be extremized is, 



r LB (/) = Ql 

v -F + 1/2 
v-F+l 



In 



F + l 



+ 



F + l/2 



v - F+l J F 
dAo ,f , d\o, v -F 



1 



dF 



dF 



C 



(27) 



where the integration constant is defined by, 

+ 1 2(i/ + 2) 



Q = jlb(F = u/2) 



u/2 + 1 v + 600/288 



-C 



(28) 

and the A derivatives arc the same as those used in go- 
ing from Equation [15] to 1171 This relaxation function 
is shown in Figure [7] and is very similar t o the corre- 
sponding function from IBarnes fe Williams! ([201 If ), their 
Figure 2. 

4. SUMMARY AND DISCUSSION 

In an attempt to better understand the evolution 
of self-gravitating collisionless systems, we have re- 
investigated a standard statistical mechanics approach 
to finding equilibria (entropy maximization), and con- 
tinued to develop a new approach (extremization of en- 
tropy p roduction), first applied t o self-gravitating sys- 
tems by IBarnes fc Williams! ([2011!) . Entropy production 
in non-equilibrium steady-state systems has been previ- 
ously investigated and found to be useful in several non- 
astrophysical systems, and so it is interesting to ask if 
the principle is relevant in gravitational systems. 

Both of our approaches use a very accurate approxi- 
mation for Inn!, and not the standard Stirling formula 
valid only for large occupation numbers. It has been re- 
cently shown that using an approximation that reflects 
correct behavior for small n, and the principle of entropy 
maximization leads to density distri butions that r esem- 
ble globular clusters (King profiles; IMadsenl 11996!) . and 
simulated pure dark mat ter halos ([Williams fc Hiorthl 
120101: [Williams etd\\2Ql($ \. Both types of systems have 
finite mass and energy, in compliance with the entropy 
maximization constraints. 

4.1. Results of Entropy Maximization 

Entropy maximization using the LB and MB statis- 
tics produces the density profiles and logarithmic profiles 



slopes, a = — dlogp/dlogr depicted in Figures S] and [5] 
respectively. 

We argue that the LB case accurately represents colli- 
sionless systems because co-habitation of the phase-space 
elements is prohibited. However, the density profiles for 
any value of v is very different from the two cosmologi- 
cal m odels, Na yarro-Frenk- White (N EW) ([Navarro et al\ 
11996! ) and N04 ([Navarro et aZ.ll2004f ). which are fits to the 
results of cosmological A-body simulations. In contrast 
to the latter, the LB density profiles have a flat central 
density slope, while at larger radii a ha s a steeper rise 
than even that of the Plummer profile ([ Schuster] 0.8831 : 
iPlummeil 119111 : iBinnev fc Tremaind 11987!) which is an 
example of a polytropic system. In Figure [5] we com- 
pare the profiles to those derived in Section IA.11 which 
are nearly identical to King models; LB profiles are poor 
matches to King models. (The vertical lines in each panel 
mark the radial position where a = 2.) Both the MB and 
LB density profiles are als o different from those found in 
iHiorth fc W illiams (|2010!) , as those models can have cen- 
tral density cusps. 

We conclude that applying entropy maximization to 
the LB case, coupled with a very accurate approxima- 
tion for Inn! produces (isotropic) density profiles that 
are distinct from any other first-principles function, or 
fits to the results of N-body simulations. Wc will return 
to this point in Section [4~3l 

The MB case is intended to represent collisional sys- 
tems because the co-habitation of the phase-space ele- 
ments is allowed. The density profiles with v = 100 and 
1000 are given in Figure [JJk, The r max scaling distance 
corresponds to the radius where a non-moving particle 
has an energy for which the distribution function dis- 
appears. The second panel of that figure ^p) shows 
the slope of the logarithmic density profile a along with 
curves corresponding to three other well-known analyt- 
ical density profiles. These MB solutions are basically 
identical to King models (see Fig ure ITU1). This result 
comes as no surprise m the light of IMadsenl (|1996j) . who 
did not use any approximation for In a;!, and obtained 
density profiles very similar to King's. In effect, aban- 
doning the Stirling approximation and treating the low 
occupation number limit with the respect it apparently 
deserves, demonstrates that the King distribution func- 
tion is the result of maximizing entropy in MB statistics, 
under the conditions of fixed mass and energy. 

The consequence of using a very good smooth approx- 
imation to In a;!, as done here, instead of the exact ex- 
pr ession which giv es discrete step-like values, as was done 
in IMadsenl (|1996f) , is that the former results in spot-on 
matches to the King profiles, while the latter display 
modest differences from the King pr ofile shape, espe cially 
for low ^(0)/cr values; see Figure 1 IMadsenl (fl996l) . 

The value of v appears to play a role analogous to 
the King scaling factor ^(0)/a 2 , where \&(0) is a rela- 
tive potential en ergy at the center, and a is an energy 
scaling constant ([Binnev fc Tremainelll987i) . Increasing 
v — or ^>(0)/a 2 — increases the concentration of the sys- 
tem. Recall that the higher the concentration, the closer 
the King model is to an isothermal sphere. 

We find that the concentrations of the King models 
that most closely resemble the Maxwell-Boltzmann den- 
sity profiles increase as the value of v increases. Recall 
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that v is the number of fine-grained micro-cells that oc- 
cupy a coarse-grained macro-cell, so as v increases the 
distribution function becomes "grainier" . At the same 
time, increasing the concentration of a King model makes 
it more closely resemble an isothermal sphere. Putting 
these two observations together implies that increasing 
the coarseness of the distribution function (by increasing 
the number of fine-grained cells contained in a coarse- 
grained cell) should result in densities that have more 
isothermal aspects. 

Finally, we note that both of the families of profiles ob- 
tained here have flat density cores, while cosmologically 
simulated halos have density cusps. 

4.2. Results of Entropy Production Extremization 

Some aspect of simulations may delay, or even dis- 
allow, the maximizing of entropy necessary to achieve 
thermal equilibrium. In such a frustrated case, it is pos- 
sible that a collisionlcss system finds a stationary state 
by extremizing its entropy production. As a more con- 
crete example of this situation, imagine a metal bar with 
one end exposed to a blowtorch and the other end in 
an ice bath. When the bar reaches a stationary state 
with a time-independent, but spatially varying, temper- 
ature distri bution, entropy product ion in the bar is an 
cxtremum (|de Groot fe MazurlH98l Ch. 5, Sec. 3). Wc 
have begun to explore the possibility that the mechan- 
ical equilibria of simulated collisionless systems can be 
explained as states of extreme entropy production. 

To this end, we have derived expressions for entropy 
production in collisionlcss systems using the Lyndcn- 
Bell statistics, and in collisional systems, using Maxwell- 
Boltzmann statistics. Further, we have found the form 
of the relaxation function required to guarantee an en- 
tropy production extreme. Recall that the relaxation 
function is the right-hand side of the Boltzmann equa- 
tion, Equation [6j that describes the coarse-grained dis- 
tribution function evolution. 

As in our previous work which was valid only for 
large phase-space occupation numbers, the Maxwell- 
Boltzmann relaxation function (pictured in Figure [6]) is 
positive for all values of F (as long as N > 1) and be- 
comes more constant as the number of particles N in- 
creases. It is clear that the basic form of the relaxation 
function is not dependent on the specifics of the distri- 
bution functio n. A visual comparison be tween Figure [6] 
and Figure 1 in lBarnes~fc William! (|2011[ ) reinforces this 
point. 

The general form of the Lyndcn-Bell relaxation func- 
tion shown in Figure [7] has similarities and differences 
with its large occ u pation number version derived in 
iBarnes fc~ Williams (|2011l) . The similarities arc that 
both the functions increase with f/r) up to just short of 
f/r) = 1, then spike towards large positive values of Tlb, 
go through a singularity at f /t)<1, then become negative, 
and eventually asymptote to Tlb = as f/r) —> 1 (see 
the bottom panel of Figure . The difference is that 
the large occupation number Tlb is independent of u, 
while its arbitrary occupation number analogue is not. 
In the latter case the values of C and v that produce the 
singularity are linked through the relationship, 



log 



100 



0.4348(C- 4.1518), 



(29) 



where v m i n is the minimum value of v that will cause a 
singularity for a given C. This expression is valid as long 
as f m in ^10- The importance of the singularity in both 
cases is that it signals that the relaxation function will 
reach zero when the coarse-grained distribution function 
/ reaches the fine-grained value 77, bringing the evolution 
of the distribution function to a halt. 

To sum up, extremizing entropy production in the LB 
case leads to a relaxation function Tlb vs. f/r) shape 
that drives T^b, and hence the entropy production rate 
<7Lb, to zero, which in turn means that the endpoint of 
evolution has a maximum (more correctly, an extremum) 
entropy. In other words, the entropy production extrem- 
ization procedure that we followed in Section 13.21 tell us 
that the final state of a self-gravitating collisionlcss sys- 
tem is a state of maximum entropy. Since in Section 13.11 
we derived such a state, using entropy maximization, it 
must be the same state. 

The important point in both the LB and the MB case is 
that extremizing entropy production leads to relaxation 
functions that drive a coarse-grained distribution func- 
tion to behave like a fine-grained distribution function. 
This has a bearing on the 'incompleteness of relaxation', 
sometimes alluded to when describing stationary-state 
collisionlcss systems. Our results suggest that if incom- 
plete relaxation does happen, it is not due to a system 
reaching an entropy production extreme. Rather, it ap- 
pears that coarse-grained evolution will proceed until en- 
tropy production ceases, when T = 0, and so full relax- 
ation will be achieved in self-gravitating systems. As an 
immediate consequence we conclude that entropy maxi- 
mization is the correct — and more direct — procedure to 
take to arrive at the description of steady-state self- 
gravitating systems. 

To test the hypothesis that for self-gravitating systems 
the state of entropy production extreme coincides with 
the state of maximum entropy, we are undertaking a com- 
parison between these analytical descriptions of entropy 
behavior in collisionlcss systems and results of TV-body 
and semi-analytical simu lations (e.g., the extended sec- 
onda r y infall model ESIM Willia ms. Babul, fc Dalcantonl 
120041 : 1 Austin et all [2001 . This further work will help 
settle exactly what role entropy production plays in de- 
termining collisionless equilibria. 

4.3. Outstanding Questions 

Aside from the future work described above, to quan- 
tify entropy production rate in simulated numerical sys- 
tems, there are also a few questions that remain unan- 
swered. 

4.3.1. Maximizing Entropy in Phase-space vs. Energy-space 

IHiorth fc Williams! (|2010f ) and the present paper max- 
imized entropy, but in different state-spaces, energy and 
(x, v) phase-space, respectively. The results of the two 
studies, for example in terms of the density profiles, 
are very different from each other. It is not surprising 
that they are different, but it is not immediately obvi- 
ous which one of the two approaches would produce a 
better description of the results of collisionlcss N-body 
simulatio ns. After looking at all the density profiles it is 
seen that IHiorth fc Williams! (120101) profiles are similar 
to those of simulated systems ((Williams fc Hiorthll2010l ; 
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iWilliams et aZ.I I2010D . while those from the LB case m 
the present work (Section 13. 11 Figures [4] and [5]) are not. 

We note that our LB e ntropy maximization results are 
what iLvnden-Behl ()1967l ) would have obtained had he not 
used the Stirling approximation. The procedure does not 
rely on any approximations, and is well motivated from 
first principles, assuming that all states are equally acces- 
sible and efficient mixing in phase-space can be achieved. 
This begs the question, do such systems exist? More 
work with simulations is needed to address this question. 

4.3.2. The Necessity of the Low Occupation Number Regime 

Our MB entropy maximization re sults of S ection IA.1I 
(and the very similar results of iMadsenl 11996) de- 
scribe some globular clusters well (lElson et ali 19871 : 
iMevlan fc Heggidll997l: IWilliams et aLll2011h . and the re- 
sults of lrliorth fc Williams! ( 2010f ) describe N-body and 



ESIM halos. These facts lead us to the following chain 
of reason. First, globular clusters, along with N-body 
and ESIM halos, are self-gravitating systems. Second, 
the entropy maximization procedures that produce suc- 
cessful models of these systems require correct treatment 
of low phase-space occupation numbers. Therefore, it 
appears that the low occupation number regime is an 
important distinction between self-gravitating and non- 
sclf-gravitating systems. 

Why is this the case? We speculate that self- 
gravitating systems are special not only because they are 
finite in extent due to the long-range nature of the force, 
but also because of correlations between spatial regions 
and energies of particles. Specifically, the center of the 
potential contains phase-space elements with the most 
bound energies, while the outskirts are populated by ele- 
ments with the least bound, or even zero energies. These 
requirements are peculiar to self-gravitating systems, and 
apparently necessitate low occupation numbers because 
the numbers of particles in these regions must approach 
zero. 

However, the above hypothesis does not explain why 
the discrete version of ln.x! must be smoothed to repre- 
sent real systems. This smoothing, though apparently 
necessary, is still not well justified. 

4.3.3. Should Maximum Entropy Imply Constant 
Temperature ? 



The steady-state of self-gravitating systems can ap- 
parently be calculated as the maximum entropy state. 
Though these systems are in mechanical equilibrium, one 
could argue that they are not in conventional thermal 
equilibrium because their kinetic temperature is mani- 
festly different across the system. Can we reconcile max- 
imum entropy defining thermodynamic equilibrium with 
the resulting non-constant temperature? (Note that a 
realistic self-gravitating system cannot have a constant 
kinetic temperature, because there is no finite-mass so- 
lution to the Jeans equation with constant a.) This con- 
tradiction is part of the reason why it made sense to 
conclude, based on ILynden-Belll (|1967l ) work that there 
is no maximum entropy state for self-gravitating systems. 

We propose that for gravitationally-bound collisionless 
systems, one must carefully separate kinetic temperature 
Tj~ from thermodynamic temperature T t . To be clear, the 
kinetic temperature we are referring to is related to the 
rms value of the peculiar velocity in a system, Tu oc (v^). 
The thermodynamic temperature is linked to the energy 
scale P that serves as a Lagrange multiplier through j3 = 
1/ {ksTt). By definition, T* must be a constant , and there 
is no contradiction as in ILynden-Belll (|1967t ) — systems 
in the maximum entropy state are characterized by a 
constant thermodynamic temperature, even as they have 
a varying kinetic temperature. 

We can also address a related question involving the 
temperature and energy change that appear in the ther- 
modynamic definition of entropy, dS = 5Q/T. This tem- 
perature must be a constant, so we would associate this 
with the thermodynamic temperature T t . SQ (normally 
called heat) refers to all non-work exchanges of energy, 
and so must be replaced by Se because in collisionless 
self-gravitating systems what is being transferred is total 
energy, potential and kinetic, not just kinetic. Entropy 
changes can occur due to changes in mass distribution as 
well as kinetic energy transport. 
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APPENDIX 
MAXWELL-BOLTZMANN STATISTICS 
Entropy Maximization 

Using the results of § 11.31 the multiplicity function for a system obeying MB statistics combined with the definition 
of entropy (Equation [1]) produces, 



■5'] 



MB 



lniV! - ^hm 4 ! + N\nv 



(Al) 



where the summation involving the macro-cell occupation rii runs over the number of macro-cells. We will assume 
that N 1, so that IniV! = N I n N — N. However, we will not use the Stirling approximation for the second term on 
the right-hand side of Equation IA11 With our approximation, Equation ll2l we can now rewrite the entropy as, 



<5mb — S] 



MB.O 



-h 



[(m + 1/2) In (ni + l) + Xo, ni ] 



(A2) 



where S'mb.o = Nk& [In (Nv) — M In 2ir/2N] is a constant, and M is the total number of macro-cells. 
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We now replace the discrete macro-cell occupation number rij with the coarse-grained distribution function / to 
produce 



Smb 



a k B 

JMB.O 

VW 



vf 



"f 



dx dv. 



(A3) 



To maximize this entropy function, we set SSmb 
that results is, 

(vf/i) + 1/2) f]d\ ^f/ v 

V 



if + 1,2 ) { if + 1 ) + A °'""" 

0, subject to constant mass and energy constraints. The expression 

fte = 0, (A4) 



(vf/r, + i) v df 



where p and ft arc undetermined multipliers associated with mass and energy conservation, respectively. The e term 
is the specific energy of a phase element located at position x with velocity v, e = u 2 /2 + $. Note that if the standard 
Stirling approximation had been employed, the second and third terms would be absent and the logarithmic term 
would simply read \n(vf /n). In this case, the usual, physically inconsistent, MB distribution function would result. 

As before, we make a change of variables, F = vf /rj, and the derivative of the A function is given by Equations 1161 
Equation IA4I can then be re-cast as, 



In (F + l) 



(264F + 276) 



576(F + 1)(F + 600/576) 



+ p + fte = 0. 



(A5) 



Again, we have not searched for an analytical solution, but graphical solutions for F for a range of e values can be 
combined to produce a plot of F(e). Specifically, Figure [8] illustrates the behavior of the normalized coarse-grained 
distribution function f jr\ for a particular value of v. Values of the Lagrange multipliers p and ft were adjusted so that 
e is always positive. Because MB statistics have no exclusion principle, f /rj can have values above 1. 



Entropy Production Extremization 
We begin by writing entropy in terms of entropy density, 



S 



ps dec, 



(A6) 



where p is mass density, s is the specific entropy, and the integral is taken over the spatial extent of the system. From 
the entropy form given in Equation IA31 the entropy density can now be written as, 



PSMB 



vw 



(F + l/2)ln(F+l) + A ,F-i ;l 



<Smb,o 



dv, 



where, as earlier, F = vf/r]. 
Taking a partial time derivative of Equation IA7I results in, 



d 



^tApsmb) = -— / ^7 



dt 



vw 



dF 



dt 



ln(F+ 1) + 



OF + 1/2) 9A , F 



(F+l) 



dF 



-B 



dv, 



(A7) 



(A8) 



where B = Smb,o/N1<;b is a constant. 

Substituting dF /dt from the Boltzmann equation into Equation IA8I results in a lengthy expression, similar to 
Equation [22] in the Lyndcn-Bcll case. We will deal with this expression term by term. For reference, the expression 
after substitution is, 



-■^ l{-vVF) 



vw 
fc B 
vw 

fcB 

vw 



In (F + l) 



(-a- V V F) 
7 In (F + l) 



In (F - 



(F+l/2) | d\ Q . 

(F + l) 
(F + l/2) 



' (F + l) 
1/2) , €>Aq, F 
dF 



dF 
dF 



B 



B 



dv 



dv 



(F+l) 



B 



dv, 



(A9) 



where 7 = {vT)/rj. 

The process to evaluate these integrals is very similar to what has been discusse d in the Lynden-Bell case. Using 
the fact that v ■ VF = V • Fv, the first integral on the right-hand side of Equation IA9I can be re- written as, 



[(F + 1/2) In (F + 1) + A ,F - FB] v dv. 



(A10) 



Note that this expression is reminiscent of the form of the entropy density in Equation IA7I 

We next turn our attention to the terms involving acceleration in Equation IA9I The fact that a is velocity- 
independent implies that a • V„/ = V„ • a/, a fact that will be used often in dealing with these terms. Let us start 
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with the first part of the the second term on the right-hand side of Equation IA91 



J (-a ■ V V F) In (F + 1) dv = - j (V„ • aF) In (F + 1) dv. 



(All) 



Integrating by parts produces two terms, one with the form, 

V„ • [aF In (F + l)] d-u. (A12) 

which equals zero after using the divergence theorem and the fact that any physical distribution function must vanish 
for large velocities. The other term resulting from the integration by parts is, 



A single component of this integral will appear as, 

l*>l,m.. g 



a ■ J (V V F) -^-j dv = a- J V v [(F + 1) - In (F + 1)] dv. 

ill appear as, 

[(F + 1) - In (F + 1)] d Vl = (F + 1) - In (F + 



dvi 



(A13) 



(A14) 



where V\ >rn i n — —Vx imax , representing the maximum speed possible for the system. Assuming that F is even in v\ ima x 
and that F(vi imax ) = so that the distribution function disappears at the maximum speed, this integration results 
in zero. Similar manipulations can be applied to the rest of the acceleration-dependent parts of the second term on 
the right-hand side of Equation I A91 resulting in the entire second term being equal to zero. Equation I A9I now has the 
form. 



dt 



(psmb) : 



VVD 



V- J [(F + 1/2) In (F + 1) + A , F - FB]v dv+ 



7 (F) 



m(F + i ) + ^ + d -i^-B 



F + l 



dF 



dv. 



(A15) 



Assuming the velocity field to be composed of mean and peculiar components v = Vo + v p , we can re-cast Equa- 
tion as, 



d_ 

dt 



(rsmb) = - V • psmbVo -—{ V • / [(F + 1/2) In (F + 1) + X ,f - FB] v p dv 

VXD I J 



7 (F) 



ln(F + l) + ^±lZ? + ^-fl 



F+l 



dF 



dv 



(A16) 



We now equate the terms in Equation lA16l to those in the continuous version of the entropy density evolution equation, 
Equation [231 

— (pa) = - V • (S + psv ) + o. ( A17) 

at 

The entropy flux S is given by the integral in the second term on the right-hand side of Eq uat ion I A 1 6 1 and represents 
randomly fluxed entropy. The remaining term is the entropy production for the system, 



vw J 



ln(F + l) + ^ + ^-F 



F+l 



dF 



dv. 



(A18) 



This equation explicitly demonstrates how the non-collisionlcss nature of the coarse-grained distribution function leads 
to changes in entropy. 

As mentioned in §[U thermodynamic non-equilibrium systems can have steady-states described by extrema of entropy 
production. We then set 5omb = 0. Taking the variation of Equation I A181 we find 



5a 



MB 



vw 



SF 



d') 
dF 



s F+l/2 dAn p 
In (F + 1) + J '- + - B 



F + l 



dF 



+ 7 



F + l/2 d 2 X 0tF 



F+l (F + l) 5 



dF 2 



dv. 



(A19) 



Since SF is arbitrary, the variation disappears only when the term in curly braces is zero. The condition for an 
cxtremum in entropy production is, 



din 7 
dF 



In (F + l) 



F + l/2 dA n ,F 



B 



F + 1/2 d 2 A 0l 



F+l (F + l) 2 



F+l dF 

The solution for the Maxwcll-Boltzmann relaxation function is, 

F+l/2 dA 0jF <5mb,o 



dF 2 



0. 



7MB (F) = Pj 



In (F + 1) 



F+l 



dF 



(A20) 



(A21) 
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where the constant can be expressed as 
P = Jmb(F = v) 



ln(j/ + 1) + 



u+1/2 



v + 600/576 



{SMB,o/Nk B ) 



(A22) 



The 7mb(-F = v) term is the relaxation function value when the coarse-grained distribution function is equal to the 
constant fine-grained distribution function value (/ = rj). 
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The set-up: statistical representations 
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Attaining final state: maximizing entropy 
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Attaining final state: extremizing entropy production 
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Fig. 1.— 



A schematic summary of the various statistical mechanical approaches to self-gravitating systems. See Section fl.4l for explanation. 
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Fig. 2. — Comparing the Stirling, IHiorth fc W illiams 1 20Tof) and Equation [9] approximations to In a;! as functions of x. The function 
ip(x + 1) = dlnxl/dx is plotted. Note the divergent behav ior of the Stirling approximation as x — > 0. (a) The raw ip functions for the three 
approximations are shown. ThepTiorth & Williams (2010) line is ver y nearly covered by the E quation |9] line, (b) Differences between the 
various approximations and Equation [9] The differences between the Hiorth & Williams (2010 ) and Equation [9] approxim ations are much 
smaller than those between either approximation and the Stirling approximation. Also, the IHiorth &c Williams! 1 120101 ) and Equation [9] 
values coincide for x = 0. 




Fig. 3. — The Lyndcn-Bcll distribution function resulting from the non-Stirling approximation. The Lyndcn-Bcll function derived using 
the Stirling approximation is also shown. 



15 
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log r 7 



Fig. 4. — (a) Logarithmic density and (b) the corresponding a = — dlogp/dlogr versus logarithmic scaled radius for the Lynden-Bell 
case. The thick solid black lines correspond to models with v = 100 while the thick dashed black lines represent models with v = 1000. 
The v ertical lines mark where the logarithmic slope of the v = 100 model is isothermal (a = 2). The a profiles for NFW. [Navarro et al.\ 
1120041 ), and Plummer models are also shown. The profiles for the v = 1000 model have been hor izontally shifted so that the location where 
a = 2 coincides with the other models. The central cusp of the NFW and lNavarro et al.\ 12004T) models (a -> 1 as r — > 0) differs markedly 
from the core present in the Lynden-Bell case (a — > as r — > 0). The presence of the core in the Lynden-Bell case is independent of v. 




Fig. 5. — Comparisons between Lynden-Bell (solid) and Ki ng ( dashed) density profiles for different v values; (a) v = 100 and (b) v = 1000. 
The King scale factors correspond to those used in Figure HOI for the same v values. The shapes of the curves cannot be brought into 
agreement for any scale factor values. 
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Fig. 6. — The behavior of the relaxation function Tmb a function of the coarse-grained distribution function / in a system that obeys 
Maxwcll-Boltzmann statistics. The various lines represent functions defined with different N values, where TV is the number of phase-space 
elements in a system. As with the LB case, as N increases, the relaxation function becomes more constant. Note that the horizontal axis 
scale differs from that in Figure [7K since there is no restriction on the value of the coarse-grained distribution function in the MB case. All 
curves correspond to cases where v = 100. 
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Fig. 7. — The behavior of the relaxation function Tlb as a function of the coarse-grained distribution function / in a system that 
obeys Lynden-Bell statistics. Panel a includes various lines representing functions defined with different C values, where C = In TV — 1 + 
M / N(u\nu — ln27r). N is the number of phase-space elements in a system, and increasing N results in increasing C. As N increases, the 
relaxation function becomes more constant. All curves correspond to cases where v = 100. As v increases, these curves develop singularities 
near / = 7], as shown in panel b. 
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Fig. 8. — The natural logarithm of the Maxwell-Boltzmann distribution function resulting from the non-Stirling approximation. The 
Maxwcll-Boltzmann function derived using the Stirling approximation is also shown. The non-Stirling approximation distribution function 
reaches zero for a finite e, unlike the exponential function derived using the Stirling approximation. 



20 



Barnes & Williams 




Fig. 9. — (a) Logarithmic density and (b) the corresponding a = — dlogp/dlogr versus logarithmic scaled radius for the Maxwell- 
Boltzmann case. The thick solid black lines represent models with v = 100 while the thick dashed black lines correspond to models with 
v = 1000. The vert ical lines mark where the logarithmic slope of the v = 100 model is isothermal (a = 2). The a profiles for NFW, 
INavarro et al. (2001), and Plummcr models are also show n. As in Figure [ 4] the v = 1000 curve has been shifted horizontally to align the 
locations where a = 2. The central cusp of the NFW and Navarro et al. (2001) models (a — > 1 as r — > 0) differs markedly from the core 
present in the Maxwell-Boltzmann case (a — > as r — > 0). The presence of the core in the Maxwcll-Boltzmann case is independent of v. 
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Fig. 10. — Comparisons between Maxwell-Boltzmann (solid) and King (dashed) density profiles for different v values; (a) v = 100, (b) 
v = 1000, (c) v = 5000. The King profile scale factor \f(0)/o" 2 is given in each panel. The dashed lines represent King density profiles 
while the solid lines show the Maxwell-Boltzmann distributions. The comparison between these two profile types reveals more similarity 
than when the[Madscn ( 199tJ) and King models are compared. 



